Simple Model of Capillary Condensation in porous media. 
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Abstract 



We employ a simple model to describe the phase behavior of 4 He and Ar in 
a hypothetical porous material consisting of a regular array of infinitely long, 
solid, parallel cylinders. We find that high porosity geometries exhibit two 
transitions: from vapor to film and from film to capillary condensed liquid. 
At low porosity, the film is replaced by a "necking" configuration, and for a 
range of intermediate porosity there are three transitions: from vapor to film, 
from film to necking and from necking to a capillary condensed phase. 



I. INTRODUCTION 

The physics of adsorption in porous media has long been appreciated for its considerable 
importance in diverse applications, including gas storage, separations and cryopumping. In 
addition, the geometry provokes intriguing fundamental questions about the properties of 
phases of the adsorbate, for which the surface curvature energy plays a prominent role |I]-[7| . 
An ongoing question pertains to the validity of numerous simple quasi-two or one dimensional 
models which omit the heterogeneity present in most such media. Moreover, there 
remains a fundamental question about the nature and existence of genuine phase transitions 
in such disordered media. 

In this paper, we attempt to answer qualitative questions about the nature of such adsorp- 
tion by employing a simple model of the geometry that we study with a simple calculational 
procedure. The geometry is shown in figure 1 and described in detail below. The procedure 
is called a "simple model" in a number of papers published over the last decade flOHTTjl. 
In the model, the system's energy is taken as a straightforward sum of bulk energy, surface 
energy and gas-solid interaction energy terms. The model is certainly oversimplified; yet 
it has provided answers to questions about similar adsorption problems which agree rather 
well with results obtained in more reliable (even "exact" ) studies of the same problem. The 
rationale for applying the present approach is that the results may yield broad trends which 
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are robust, at least qualitatively. In the present case, our results are phase diagrams for 
which a key parameter is the porosity $ of the system. As expected, the physics of pore 
filling involves a competition between adhesive and cohesive forces. As such, the results of 
our calculations are sensitive to assumptions we make about the gas-surface interaction. In 
the present case, both these assumptions and the calculational results are analogous to those 
found in the problem of wetting transitions. |[ 



II. THE MODEL 

We propose a simple model consisting of infinitely long, solid, parallel cylinders of radius 
R = 30A. The centers of the cylinders form a square lattice and the distance between sites, 
5*, is related to $ by the equation 

The configuration is depicted in Fig. [1], where is also shown the triangular unit cell consid- 
ered in the calculations. The interaction between the adsorbed atoms and the substrate is 
approximated by the sum of the potentials from the four cylinders that are closest to the 
unit cell. The potential due to each single cylinder, V, was constructed assuming cylindrical 
symmetry so that it is a function of the distance to the axis only (r). The following function 
V(r) reproduces both the potential of a flat substrate (which has often been modeled as 
4C|/ (27D 2 z 9 ) — Cs/z 3 , z being the distance from the surface) for r — > R and the potential 
of an infinite wire for r ^> R, 

V(r) = C - (2) 

K> 21D 2 {r-Rf {C/C z -R 2 + r 2 ){r-Rf' K) 

where D and C3 are the well depth and van der Waals parameters respectively |19|J20[1 and 



C = 9^CsR 2 . We adopted values displayed in table | that are intermediate between the 
strongly attractive graphite and the weakly attractive alkali metals. The form of V(r) and 
the total interaction considering four cylinders are plotted in Fig. |2|. 

Atoms adsorbed in this material below saturated vapor pressure are expected to form 
either a film adsorbed on the wall of the substrate, or a condensed phase (C) filling all the 
free space. Which of the two phases is stable can be determined by evaluating and comparing 
the free energies. If the film is stable, and the porosity is moderately low (<50%), the fluid 
may form bridges or "necks" between neighboring cylinders (see Fig. |3|). We will refer to 
this configuration as "necking" (N), and apply the term "film" (F) to the case where no 
necks are formed. The letter E refers to the empty pore configuration. This actually means 
the presence of a low density vapor phase, which is "empty" only by contrast to the higher 
density liquid phase. 

For a given value of the chemical potential \x below saturation (/io), we evaluate the 
grand free energy per unit length 

Q = F/L - fiN/L. (3) 

F is the Helmholtz free energy, consisting of the substrate-fluid interaction, a bulk free energy 
and the surface energy. Assuming translational symmetry along the coordinate parallel to 
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the axis of the cylinders, z, the grand free energy per unit length of the fluid contained in 
the unit cell is, for the phase C 



tt/4 f S/(2cosfo)) 
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Jr. 



n c = n I d<p / V(r)rdr - n(ji - fM )(S - rtr min )/8 + cr-r min . (4) 
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Here a and n are the bulk surface tension and number density of the fluid at a given 
temperature, and r min is a cutoff distance, chosen as the radius where the potential due to a 
single cylinder is minimum. We are taking the surface tension for the solid-liquid interface 
to consist of the sum of the liquid-vapor interfacial tension and the integrated solid-fluid 
potential energy. In this model, we are assuming that a fluid of density n fills all the 
available space homogeneously , except the region r < r min . For the phase F we consider a 
homogeneous film of density n formed between r m i n and r e (<p), the energy reads 



\n{n - // ) fo /4 d<p{r e (<p) 2 - r 2 min ) + o\ r min , 



where r e (<p) is the equilibrium profile, and 4>n is the angle subtended by the neck (see Fig. 

|3p if present, taken to be zero for a film. A prime in this expression means a derivative with 
respect to (p. 

By minimizing the energy Qp, we get a differential equation for r e {4>): 

^-^-V[r e (<f ) )} = ^ M y (6) 

where R[r e (4>)] is the radius of curvature, 

R[rU)\ = Wl!±lWg (7) 

Notice that if the surface tension is small the right hand side of eq. ^| can be neglected 
and the following condition results 

V[r e ((f>)] = n - mo- (8) 

This is the equation for the equipotential curves; its use corresponds to the venerable Polanyi 
theory of adsorption pi| . In the present study, we use this relation, instead of solving 



eq.|, because it greatly simplifies the calculation. The approximation is good in a range of 
chemical potential such that 

fio - ii >| a/nR{r e ((/))} | (9) 

for all 0. It will be shown below that in the case of He, the relevant values of /i, where we 
find phase transitions, fulfill this condition. For Ar, that is not the case at temperatures 
close to the triple point (i.e., when the surface tension is maximum). These behaviors will 
be discussed in the next section. 
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III. RESULTS 



In Fig. [| a) we display the grand free energy of the phases F and C for Ar at T = 130/^ 
and porosity $ = 0.42, which corresponds to a separation d = 10A between the cylinders. 
Observe that Q c an d &f > for 5li* = (/i — /io)/ e < — 1-7; hence the empty phase is 
stable in this regime. However Qp < &>c and Qp < for —1.7 < 5fi* < —0.28, so the film 
is stable there, while for 8li* > —0.28 the completely filled phase is stable. In Fig. |] b) 
we plot the corresponding isotherm, together with the derivative dN/d/i. The steps in the 
isotherm correspond to the transitions E— >F and N— >C, and the peak in dN/dii corresponds 
to the transition F— >N . From similar analysis for different porosities, we construct the phase 
diagram shown in Fig. [5[ The line of the phase diagram is dotted in the region where the 
condition given by eq. |9] fails, according to the following criterion: in the case of films the 
quantity cr/nR[r e ] will be smaller than a/nr min = —bii* min , so we consider that the condition 
(0) is fulfilled if 5ii* < 10 5^^. In the case of necking, the curvature radius for (f> ~ (fi N is 
not necessarily larger than r min and therefore r eq (<p ~ <f)^) is poorly determined. However, 
the error introduced to the energy by this difference in the profile is expected to be small. 
For temperatures closer to the triple point, the surface tension (and therefore — 5/i^ in ) is 
larger, and the transitions occur at values of ii such that 8ii* > 10 5{i^ in , so the analysis 
done using eq. |S] becomes meaningless for such low temperatures. As is seen in Fig. ^ the 
phase N is possible only for very low porosities (small d), where the volume- filling "price" 
of necking is smaller. The upper part of the diagram is qualitatively equivalent to the one 
obtained for a slab with large separation between walls, (see Fig. 4 in ref. ||10|| ) 

The corresponding energies calculated for He are shown in Fig. [| for two different 
temperatures and $ = 0.42, together with the components of Qp, on the right. By means of 
the same procedure described above for Ar, we derived the phase diagrams displayed in fig. 
[7|. The scenario is the same as for Ar. Comparing the results for T = K and 3 K, we notice 
that at higher temperature there is a wider region where F is stable. This can be explained 
in terms of the values of the surface tension: for T=3 K, a/n is 30% smaller than for T=0 K. 
This is evident in Fig. || on the right, where we can see that the contribution of the surface 
energy to flp is more important for T=0 K than for 3 K. As a consequence, the phase F, that 
requires the formation of an interface, is more favorable at higher temperatures. For Ar, the 
region of the phase diagram corresponding to the phase F is even narrower, as is expected 
since the surface energy relative to the substrate potential energy (in relative terms) is larger 
for Ar than for He at any temperature. 

IV. COMMENTS 

In this paper we have derived the phase behavior for gases exposed to an array of parallel 
cylindrical strands of material. Using a number of simplifying assumptions, this calculation 
yields a kind of generic behavior, exemplified by the similarity of the behavior predicted for 
a classical fluid (Ar) and an extreme quantum fluid ( 4 He). The behavior is such that low 
porosity ($ < 0.4) geometries exhibit a necking transition at low chemical potential, followed 
by capillary condensation. For a narrow range of intermediate porosity (0.4 < $ < 0.5), the 
behavior is different: a thin film forms, followed by necking (with a singular derivative but 
no discontinuity in the isotherm), followed by capillary condensation. At porosity $ > 0.5, 
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there is a transition from film to capillary condensation and necking does not occur at 
any chemical potential. We note that the Kelvin equation does not agree with the present 
results for the threshold of capillary condensation. This discrepancy is well known from 
similar studies of small pores, due to the equation's neglect of the substrate potential and 
the resulting depleted density region near the wall (so that the pore radius is ill-defined). 

m 

This research has been funded in part by the Army Research Office, the Petroleum 
Research Fund, the American Chemical Society, CONICET and the University of Buenos 
Aires. 
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TABLES 



TABLE I. Well depth and van der Waals coefficients D and C3 used in the calculation and 
parameters of the gas-gas interaction, e and a from ref. [23|. 





e(meV) 


(7(A) 


D(meV) 


C 3 (meVA 3 ) 


Ar 


10.34 


3.40 


50 


1000 


He 


0.88 


2.56 


7 


150 
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FIGURES 

FIG. 1. Schematic diagram of the configuration; the unit cell is shown with dashed lines. 
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FIG. 2. a) Equipotential curves due to the four closest cylinders (with = 0.42), for He. 
Contours are labelled by potential energy values, in meV. The dashed line shows a boundary of 
the unit cell, b) Potential contributed by a single cylinder for He (dashed) and Ar (solid). 
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FIG. 3. Schematic views of the film phase (left) and necking phase (right). 
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FIG. 4. a) Grand free energies for Ar at T = 130-fT (<& = 0.42). b) Number of atoms per unit 
length (solid line) and dN/dfi (dashed line, in arbitrary units). The vertical line at the isotherm 
step corresponds to a delta function in this derivative. 
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FIG. 5. Phase Diagram for Ar at T = 130K. The dotted portion of the curves are explained 
in the text. 
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FIG. 6. Pore-filling behavior for helium at T=0 and 3K. Left panels show the grand free energies 
CIf (solid line) and Q,c (dashed line) for He at the temperature indicated. Right panels show the 
energy contributions to Qf, due to the surface (dotted line), substrate (dashed line) and volume 
(dotted-dashed line). 
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FIG. 7. Phase Diagram for He, at the temperatures indicated. 
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